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We study quantum spin ice in an external magnetic field applied along a (100) direction. When 
quantum spin fluctuations are weak, elementary excitations are quantum strings with monopoles at 
their ends manifested as multiple spin-wave branches in the dynamical structure factor. Strong quan- 
tum fluctuations make the string tension negative and give rise to the deconfinement of monopoles. 
We discuss our results in the light of recent neutron scattering experiments in Yb2Ti2C>7. 

PACS numbers: 75.30.-m, 75.40.Gb, 75.40.Mg 



The quest for novel quantum phases and elementary 
excitations is one of the central themes in condensed- 
matter physics. The notion of an elementary excitation is 
conventionally associated with a point-like object, as the 
term quasiparticle implies. A natural question is whether 
elementary excitations in quantum materials could re- 
semble strings, rather than particles. String excitations 
were recently found in spin ice Dy 2 Ti 2 07 [U [2], a frus- 
trated ferromagnet with fractionalized excitations known 
as magnetic monopoles [21 S]. In an applied magnetic 
field, excitations are strings of misaligned spins connect- 
ing two monopoles of opposite charge. 

Conventional spin ice is a classical magnet with Ising 
spins [5]. Therefore, magnetic monopoles and strings in 
it are classical objects whose dynamics are due to ther- 
mal fluctuations. In this letter, we propose that string 
excitations with inherent quantum dynamics may exist 
in quantum spin ice, a new family of spin-ice materi- 
als exemplified by Tb 2 Ti 2 7 and Yb 2 Ti 2 7 [6][7]. In 
these compounds, spins exhibit substantial quantum fluc- 
tuations. We demonstrate that, in a certain regime of 
coupling constants, elementary excitations of quantum 
spin ice are strings with quantum dynamics. The calcu- 
lated dynamical structure factor S(u>, k) reveals multiple 
branches of excitations that correspond, loosely speak- 
ing, to strings of different lengths. As the applied field 
increases, these branches gradually separate and the low- 
est one evolves into a magnon. We connect these findings 
to recent experiments on neutron scattering in Yb 2 Ti 2 C>7 

mm- 

We begin with a toy model of quantum spin ice on 
the two-dimensional checkerboard lattice, Fig. [T] The 
point of departure is classical spin ice, in which spins have 
projections S? = ±1/2 on local directions shown in 
Fig. [T^i. Magnetic charge on a crossed plaquette (planar 
tetrahedron) is defined as = — X^eKl w ith = 
±1 for sublattice A (B). The ground states of the classical 
spin-ice Hamiltonian, 

H = J2 J2 ^f^=E J ^/2 + const, (1) 

obey the Bernal-Fowler rule, — 0, on every tetrahe- 
dron [5], Next we apply a weak magnetic field in the ac 
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FIG. 1: (a) The checkerboard lattice. A and B denote two 
symmetrically inequivalent planar tetrahedra, and arrows the 
local Zj directions, (b) The fully-polarized state when the 
field is applied in the c direction. Arrows denote the spin 
orientations, (c) A string of flipped spins (light green) binding 
a Q = +1 monopole (red solid circle) and a Q = — 1 one (open 
blue circle), (d, e, f) -ImS aa (oj, k) for k b = 0. B/h = 0.5, 
1.5, and 4.5 respectively. 



plane. In the local frames, the perturbation reads 

H 1 = -J2(hS* + B m St). (2) 

i 

Here we chose the local y-axes to be orthogonal to the 
field and introduced cosines rji = c-Zj = ( — \) Ci j\pi. The 
Zeeman term ^ has two effects. Its longitudinal com- 
ponent B breaks the degeneracy of ice states and favors 
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FIG. 2: Definition of X c operators. The operator X2 changes 
the orientation of the string segment with c — 2 from 1 to r in 
state (a) and results in a new state (b). The operators X-i $ s 
fall outside the range of the string and act trivially on (a). 

a fully magnetized state, Fig. [T]d. The transverse compo- 
nent h induces quantum fluctuations of spins. We treat 
B and h as independent parameters in the toy model. 

Flipping a single spin in the fully magnetized state cre- 
ates two monopoles with Q = ±1, which can be pulled 
further apart. The process creates a string of spins 
aligned against the field and connecting the monopoles, 
Fig.[TJ:. For h = 0, the energy of a string with n segments 
is J+Bn/^/2. For weak fields, the Hilbert space thus sep- 



arates into near-degenerate subspaces with a fixed num- 
ber of strings. The transverse part of the Zeeman term 
([2]) mixes states in the same subspace through quantum 
tunneling, inducing quantum motion of strings. We use 
degenerate perturbation theory in the subspace with a 
single string to construct an effective theory of its quan- 
tum dynamics. 

The shape of a string is specified by its segments 
{si,S2 . . . s„}, or {s^ for short, which take on the val- 
ues r = (0,1,1) and 1 = (0,-1,1) in the a&c-frame. 
The string thus propagates upwards in Fig. [TJj from the 
Q = +1 monopole at s + to the Q = — 1 monopole at s_. 
Because of the constraint s_ — s + = Yli=i s *> the state 
of a string is fully specified by its shape and location of 
one of the ends, |s + , {s}). We introduce a hybrid basis 
with fixed shape {s^}, c-coordinate of the monopole c+, 
and the 6-component of the total momentum kf,: 

\k b ,c + ,{s l })=Y / e lkb{b++b - )/2 \b + , (3) 
K 

To the first order in h, the motion of a string involves 
removing or adding a segment at one of the ends, with 
an effective Hamiltonian 



H oS |c+, {si . . . s„}) = (J + nB/V2)\c + ,{ Sl . . . s n }) - (h/2)e lk » K/2 \c+, { Sl . . . s„_i}) - (h/2) e - ik ^ 2 \c + + 1, {s 2 . . . s„}> 
- (ft/2) ]T e-^ b "+i/2| c+; {si . . . Sn+l}) _ {h/2) J- e ik » b °/ 2 \c + _ 1, {so . . . Sn} ). (4) 

8 n +l S 



Here bi stands for the b-component of the vector Sj. We 
have omitted the momentum index to simplify the nota- 
tion. 

When kb = 0, diagonalization of H c g is simplified by 
the presence of multiple reflection symmetries. Define 
the parity operator X c that switches between the 1 and r 
orientations of the segment with coordinate c and keeps 
all other segment variables Sj intact (Fig{2j), e.g. 



3 c— c+ 



.1. 



= |c+,. 



-r ■ 



(5) 



When X c falls outside the range of the string, c + < c < 
c_ , it acts on the vacuum state, which is symmetric, so we 
set X c |c + , {s}) = +|c + , {s}) in this case. It can be seen 
that Xl = 1 and [X c ,xj\ = 0. Although X c does not 
preserve the coordinate of the other end of the string s_ , 
at ki, — its horizontal displacement makes no difference; 
therefore, [X c ,H e s] — 0. Thus, all fcb = eigenstates of 
H c ff can be classified by their parities under {X c } and 
H e g becomes block-diagonal. The most important states 
have all even parities, X c = +1. An all-even state of a 



string of length n and longitudinal momentum k c is 
\k c ,n)=2~^Yl E e 4Mc++C - )/2 | C+ ,{si...s„}>. 



c+ si...s n 



For them, the Hamiltonian Q simplifies, 



H cS \n) = I J 



^ I In) — %/2/icos — 
y/2j 2 



E 

m— n± 1 



(6) 

\m). (7) 



The above is equivalent to the one-dimensional prob- 
lem of a particle on a lattice subject to a constant force 
—Bj\j2 and a hard wall at n — 0. For B -C h, we use 
the continuum approximation to find the spectrum: 



Ej{k c ) = J-2V2 



h cos 



k c 



+A,- 



n>2 j 

— — B /icos — 
2 2 



1/3 



(8) 



Here Xj are roots of the Airy function. When B ^S> h, 
the lowest eigenstate is a single misaligned spin with the 
dispersion 



Ei(k c ) = J + 



B 



V2h 2 
B 



(1 + COSfc c ). 



(9) 
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Likewise, H e g can be diagonalized in odd-parity sectors 

Strings can be directly observed in neutron scatter- 
ing experiments. A scattered neutron flips a spin in the 
fully-polarized background, creating a string of length 1. 
The intensity of scattering is proportional to the over- 
lap between a length- 1 string and a string eigenstate of 
H c g squared. Fig. [I] shows the dynamical structure fac- 
tor — lmS aa (u>, k) at several values of B/h for kf, = 0. 
For this direction of k, the spectral weight comes solely 
from states with all-even parities, X m = +1. For B < h, 
the spectrum consists of overlapping bands, whereas for 
B h the bands separate and the spectrum becomes 
dominated by the shortest string consisting of a single 
flipped spin, in essence a magnon. 

For general k, we used the Lanczos method to calcu- 
late the spectrum numerically and found similar behav- 
ior. Parities X c are no longer good quantum numbers; 
therefore, more bands appear in the spectrum. 

The case of three-dimensional quantum spin ice, with 
5 = 1/2 spins on the pyrochlore lattice, proceeds along 
similar lines. The most general exchange Hamiltonian is 
written in local axes (Fig. [3^,) as [TU] 



(d)B/Jz 



JzzS^Sj — J z ±[Sf(£ijSj' + CijSj ) + ( 



(ij) 



-J±(S+SJ + h.c.) - J ±± (QS+S+ + h.c.) (10) 

Here Qj = Qi are phase factors, and i and j labeling 
spin sublattices to 3. Specifically, Coi = C23 = — 1, 
C02 = C13 = cxp(i7r/3), C03 = C12 = exp(-wr/3), and 
Qi = 0. The J zz term describes classical spin ice, whereas 
the three remaining terms create quantum fluctuations. 
A magnetic field applied in the [001] direction adds the 



Zeeman term 



"0,3 = —"1,2 



70,3 = —7l,2 



9xy 



1 

V3' 



A),3 = — Pl,2 = 



with the cosines 

9xy 
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Z V2' 



(11) 



where g xy and g z are the principal components of the g- 
tensor. In what follows we assume that the spin-ice term 
J zz dominates and treat the rest of the terms as pertur- 
bations. The z Zeeman term favors the fully-magnetized 
state (Fig. ^ijp) . Excitations are open strings connecting 
a pair of monopoles with Q = ±1. Magnetic charge is 
defined as usual, = ~ e Kl SieH &ii where M stands for 
a tetrahedron and eg = ±1 for tetrahedra of sublattice 
A (B). 

The state of a string |s + ,{s}) is again parametrized 
by the location of its Q — +1 end s + and by its shape 
{s} = {si, s 2 . . . s„}. String segments Sj have four pos- 
sible orientations: bo = (1,1, l)/4, bi = (— 1, l,l)/4, 
b 2 = (1,-1, l)/4, and b 3 = (-1,-1, l)/4. A segment 
with orientation bo or b3 must be followed by a segment 
with orientation bi or b2, and vice versa. 



\ * IK-AM) 
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FIG. 3: (a) A and B denote two inequivalent tetrahedra in 
the pyrochlore lattice and ~ 3 four sublattices. The gray 
and black arrows show the local x and z directions. The 
abc vectors specify the local frame for one sublattice. (b) 
The fully-polarized state when the field is applied in the c 
direction. Arrows show the spin orientations, (c) A string of 
flipped spins (light green) binding a Q — +1 monopole (red 
solid circle) and a Q = — 1 one (blue open circle) (d,e,f) The 
neutron scattering spectra for the momentum transfer k || B. 
B/J z ± — 1,3, and 6 respectively. 



The effective Hamiltonian in the subspace of a single 
string is 



-VZJ Z± K 1 - J ± K 2 - 2J ±± K 3 + V (12) 



Kinetic terms K\ and K2 describe first and second- 
neghbor hopping of the string ends, whereas K 3 describes 
the hopping of a string of length 1. V = J + for 
a string of length n. The explicit form of Ki is given in 

[? ]• 

Fig. [3] shows the neutron scattering spectrum 
— (Im5 aa + lmS bh ) calculated with the aid of Lanczos 
diagonalization, for momentum transfer k || B [? ]. 
We set J± — J±± = 0.36J Z ±, and g X y/g z — 2.4 as 
in Yb 2 Ti 2 07 [5]. The spectral features resemble those 
of 2D strings (Fig. [I]). The branches gradually sep- 
arate as the string tension increases with B. When 
J z ± oj J ± r~j J ±± <c -B <C J zz , the monopole dynam- 
ics is dominated by the x and y Zeeman terms whereas 
the string tension is provided by the z term. 
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FIG. 4: Loop-flipping processes in (a,b) checkerboard lattice 
lr) o jrl) and (c,d) pyrochlore lattice jb^babi) o jbibsb2). 



To the first order in perturbations J. 



z±,±,±±, 



trans- 



verse fluctuations induce the motion of a string's ends. 
At higher-orders in these couplings, the string's shape 
can change as well. The process involves the simuta- 
neous reversal of spins around a closed loop (minimal 
length 4 in square ice and 6 in pyrochlore ice) jTTJ . 
In square ice, a state | . . . lr . . .) turns into | . . . rl . . .) and 
vice versa, Fig. [4j When the position of the monopole 
and the anti-monopole are both fixed, these fluctuations 
can be mapped onto aS = 1/2 XY chain [13] . with spin 
values t z = ±1/2 representing r and 1 segments, and the 
Hamiltonian 



state of the 2D quantum spin ice is likely another con- 
fined phase separated from the fully-polarized state by 
the transition at B c . 

In the quantum spin-ice material Yb2Ti2C>7, the cou- 
plings associated with quantum spin fluctuations, viz. 
Jz± , J± ! and J±± , are comparable with the spin-ice term 
Jzz [S]- Therefore, perturbative calculations don't ap- 
ply to it directly. Nonetheless, the physical picture is 
expected to hold beyond the perturbative regime if the 
material lies in the phase that is adiabatically connected 
to the magnetized state. A recent experiment indicates 
the ground state of Yb2Ti2C>7 is a ferromagnet [9 . The 
spontaneous magnetization in a (100) direction acts as a 
"molecular field," creating nonzero string tension even in 
the absence of an external field. We expect that strings 
in quantum spin ice can be detected by neutrons and 
photons. It would be particularly interesting to observe 
a continuous evolution of string excitations in an increas- 
ing magnetic field applied along a (100) direction. 

The authors would like to thank Rudro Biswas and 
Martin Mourigal for useful discussions. Research was 
supported by the U.S. Department of Energy, Office of 
Basic Energy Sciences, Division of Materials Sciences and 
Engineering under Award DE-FG02-08ER46544. 



n-1 
i=l 



H.c), 



(13) 



where VjjD = 0(h 4 /J 3 ). Quantum fluctuations reduce 
tension of the string to B/y/2 — 2|V2d|/t- When the ap- 
plied field is below the critical strength B c — 2v / 2|V 2 d|/t'", 
the energy cost for string excitations is negative and the 
fully-polarized state becomes unstable. A similar tran- 
sition occurs in the pyrochlore quantum spin ice, where 
a string can be mapped onto an XY chain with second- 
neighbor interactions only [Fig. Qc) and (d)] , 



flue 



n-1 
i=l 



2i- 



H.c). (14) 



The string tension is reduced by 2|V3d|/tt- When B is 
below the critical value B c = 2v / 3|V'3d|/7'", the polarized 
state becomes unstable. 

The fate of the ground state below B c depends on the 
dimensionality. On the one hand, the zero field ground 
state of the pyrochlore spin ice in the perturbative regime 
Jz±,±,±± <C J Z z is a U(l) spin liquid with deconfincd 
monopoles|12j. Therefore, the transition at B c could 
be associated with deconfinement of monopoles. On the 
other hand, given that the compact quantum electrody- 
namics is always confined in 2D [2], the B = ground 
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Supplementary Material 



Even and odd parity states of quantum strings in 2D 

We define the following reflection operators X c (Fig. 2): 

X c |s + ,{s 1 ...s m ...s n }) = ( |s+>{bi ••*»}); m = c- C+ + l 

c| +' 1 m nSI \ |s+,{si...s m ...s n }); c< c+ orc> c+ + n- 1 v ; 

Here we have introduced shorthand notation 1 = r and r = 1, and c + is the c component of the monopole position 
vector s+. X c reflects the orientation of the segment located at c. If no segment is located at c, X c effectively acts on 
vacuum, and therefore we define that X c acts trivially on such a state. It can be seen that X% = 1 and [X c , X c >] = 0. 
Now we consider the subspace of kb = 0. The effect of X c acting on \kb = 0, c + , {si}) is given by: 

yiL_ n , /- , _ f |fc 6 =0,c +) {si...s TO ...s„}); m = c-c+ + l 

X# fc -0,c + ,{s 1 ...s m ...s„})- j | fcb=0jC+){si Sm Sn}) . c<c+0 rc> C+ +n-l ^ 
Hence, fc& = states form an invariant subspace of X c operators. When kb = 0, the effective Hamiltonian (4) becomes: 
H cS \c + , {si . . . s„}) = (J + nB/\/2)|c + , { Sl . . . s n }) - (h/2)\c + , {s x . . . s n _i}> - (V2)|c+ + 1, {s 2 . . . s n }) 

- (V2) J] l c +' { § i • • • s »+i» - (V2) £ |c+ - 1, {s . . . s„}). (S3) 

Sn+l s 

We have dropped the kb index for clarity. It can be seen that [X c , H e g] = in the kb = subspace. In what follows, 
we present eigenstates and eigenvalues of the above Hamiltonian in different parity sectors. 



All-even parity states 



We consider the all-even states for which X c = 1 for all c: 

\c+,n) =2-"/ 2 £ |c + , Sl ...s„) (S4) 

si . . .s n 

Acting H c ff on |c+,n), 

#eff|c+, n) = (J + r ^)\c + ,n) - -j=\c+,n+ 1) - -^=|c + ,n- 1) - -^|c+ - l,n + 1> - -^|c+ + 1) (S5) 

The effective Hamiltonian (|S5| is reduced to a one dimensional two-body problem where one particle is located at c + 



and the other is at c_ = c+ + n > c + (Fig |Sl| . Note that (S5) is translationally invariant in c, and the states with 
momentum k c are henceforth constructed (6): 

\k c ,n) =Y,^ (c++c - )/2 \c+,n) =2- n / 2 J2 E e ife ^++ c -)/ 2 |c + , Sl ... Sn ) (S6) 

H e g becomes block-diagonalized as shown in Eq.(7): 

H cS \k c ,n) = (J + ^|)l fc c,n) - V2hcos{^)\k Cl n + 1) - V2h cos(^)|fc c , n - 1) (S7) 
We shall not repeat the solution in the all-even sector, which has been presented in the main text. 



Odd-parity states 

Now we discuss odd-parity sectors. We consider states that are odd under one X c operator and even under all 
the others. It doesn't exhaust all possibilities, yet our discussion on this special case can be generalized to more 
complicated situations straightforwardly. 



G 




D.4 0.6 



FIG. SI: (a) The effective one-dimensional problem for the all-even sector. The string (purple rectangle) binds a monopole (red 
dot) and an anti-monopole (blue dot), (b) The effective one-dimensional problem for a odd-parity sector, (c) The dispersion 
of the first few even (blue solid line) and odd (red dashed line) strings modes. B/h = 0.5. 



Without loss of generality, we impose the condition: Xo = — 1, and X c — 1 for c ^ 0: 

\c+,n) =2-"/ 2 J2 (-) x(1 " c+) |c+,s 1 ...s„) l-n< C+ <0 (S8) 

si . . .s n 

x(m) = 1(— 1) for s m = l(r). We mark such a sector as {odd,c = 0} where c stands for the position where the 
odd-parity condition is imposed. Acting H e g on |c+,n), we yield 

H eS \c+,n) = {J + nB/V2)\c+,n) - (h/V2)( ^ \c+, m) + |c+ - 1, n + 1) + \c + + 1, n - 1)) (1 - n < c+ < 0) 

m— n±l 

(S9a) 

H cS \0,n) = (J + nB/V2)\0,n)-(h/V2)( J2 |0,m> + | - 1)) (n > 2) (S9b) 

m— n± 1 

i? e ff|l — n, n) = ( J + U.B/ \/2)|l — n, n) - (/i/v^)(|l -n,n + l) +| 1) + \2-n,n- 1)) (n > 2) (S9c) 

F eff |0, 1) = ( J + B/V2)|0, 1) - (h/V2)(\0, 2) + | - 1, 2)) (S9d) 



Once again, Hamiltonian ( S9 1 is equivalent to a one dimensional two-body problem in which one particle is located at 
c + and the other at c_ = c + + n > c + . We see that the odd parity condition at c = effectively impose a constraint 
on the particle coordinates; the monopole is constrained to c+ < whereas the anti-monopolc c_ > 1 (Fig SI). The 
monopole and the anti-monopole are localized near c = due to the binding force between them. Therefore, in terms 
of string, the odd-parity condition imposed at c effectively pins the string at the very same location. This is in contrast 
to the all-even sector where strings are free to move. 

In the limit of B <C h, the continuum approximation gives the energy spectrum 

2\ ' 



E tJ = J — 2V2h - B/V2 + (A, + \ 3 ){hB 2 )i/V2 (S10) 

where A.; are zeros of Airy function. Different from eigenstates in the all-even sector, the eigenstates in the odd 
sector {odd,c\ are labeled by two indices each corresponding to one end of the string. In the limit of B h, 
perturbation theory gives the energy of the lowest eigenstate 

E afi = J + B/V2-V2h 2 /B (Sll) 

The eigenstates in two different sectors {odd,c} and {odd,c'} are related through rigid translation. Therefore, we 
can construct the Bloch states with crystal momentum k c : 

\k c ,{i,j}) = 52e ik < r \{i,j},r) (S12) 
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where \{i,j},r) is the eigenstate {i,j} in sector {odd,r}. Since the odd-parity states are all localized, the band 
dispersion of the Bloch state is flat: 

E iJ (k c )=E iJ (S13) 
in contrast to the dispersive bands of all-even states. 

The explicit form of Ki 

We introduce the momentum basis for strings in pyrochlore quantum spin ice: 

|k,{sJ)=^ e ^( s + + s -)/ 2 |s + ,{sJ). (S14) 

Here s + and Sj have been defined in the main text. For a string with momentum k, 
JTiKsi . . . s„}> = 7 *(s n _i, s„) e tk - s »/ 2 |{ Sl . . . sn-i}) + J2 l(s n ,s n+1 )e-^ s ^/ 2 \{ Sl . . . s n+1 }) 

S n +1 

+ 7 *(s 2 , Sl )e- ik - Sl / 2 |{ S2 . . . s„}) + 7(si, So)e lk - s °/ 2 |{s . . . s„}) (S15a) 

so 

K 2 \{ Sl . . . S„}> = e *-(s„_ 1+ s„)/2| {si _ _ _ Sn _ 2}) + e -ik.(s„-s„)| {si _ _ _ §n}) + J- e -*-(Sn +1 +s n+2 )/2| {si _ _ _ Sn+2}) 

S„ + l,S„ + 2 

+ e -ik.(s 1+ s 2 )/2| {s3 ___ Sn})+e ik.(s 1 -s 1 )|{s 1 ... S7i} ) + e jk '( s - 1 + s °)/ 2 |{s_ 1 ...s n }) (S15b) 



S-l,So 



^ 3 |b , 3 ) = £ C [ o,3]^os(k. b °- 3 + bt )|b I ); K 3 \b h2 ) = Cn.2i.cos (k ■ bl2 + b M |b,). (S15c) 



s=l,2 x 7 i=0,3 



We have omitted momentum index k for clarity. Here 7(bj, bj) = jij is a 4 x 4 anti-Hcrmitian matrix with nonzero 
elements 701 = i - re""/ 3 , 702 = e* 57r / 6 - re-' l7r / 3 , 713 = e"/ 6 + re"/ 3 , 723 = i + re l7T / 3 ; and r = g xy B/('sV2g z J z ±). 
The short- hand notation s, means exchanging b -H- b 3 and bi «-> b 2 . 

Calculation of the dynamical structure factors 

In this section, we describe the technical details of calculating the dynamical structure factors ^""(k, ui), a = a, b, c. 
The imaginary part of the dynamical structure factor is given by Kubo formula, 

- lmS aa (k, W )/tt = £ |/£(k)| 2 <5(w - £?„ + £7 G ). (S16) 

n 

The summation is over all excited states n with energy E n . The scattering amplitude /"(k) is given by 

/«(k) = Ee* R (n|^|G) (S17) 

R 

Here |G) is the ground state. R is the spatial position of operator S^. 

We discuss the 2D quantum spin ice at first. Here we just consider S aa for simplicity. Within the leading order in 
h/J, the ground state is the fully-polarized state. The effect of on the ground state is to flip the spin at R and 
equivalently to create a string of length 1. Therefore, 

$>£e* k ' R |G) = ^(e lk - 1/2 |k,l) + e* r / 2 |k,r))^ |k» (S18) 

R ^ 2 

up to a normalization constant. Here |k, 1/r) is the length 1 string state with total momentum k and bond orientation 
1/r. The low energy scattering is dominated by the scattering processes between the ground state and single string 



states. Therefore, we sum over all single-string states instead of the whole Hilbert space. The scattering amplitude is 
determined by the overlap between |k,-0) and the single-string eigenstate |k, n) of H c ff\ 

/„(k) = <k,ri|k,V> (S19) 
We obtain the formula for calculating —ImS aa (k,uj) in practice: 

- -ImS aa (k,u) = --Im(k,V| -— ^ „ |k,V)- (S20) 

7T 7T UJ + lU + — H c g 



The standard Lanczos technique is used to evaluate (S20|. Note that H B s is block-diagonalized in k and only the 
states with momentum k are needed. For each k, A max = 100 ~ 150 Lanczos iterations are used to construct Lanczos 
basis vectors, starting from the initial state |k, ip) . The string states are truncated at maximum length L max — 15 ~ 20 
due to the limitation of computer memory size. The convergence is checked by varyi ng A max and L max . A finite but 
small broadening factor r is used to replace the infinitesimal positive number + in (S20). 

For the special case kb — 0, the calculation of (S20) is greatly simplified. The initial state |k, ip) becomes 

|k,V) - ^(IM + l k , r » = \k c ,n= l) cvcn (S21) 



which is an all-even state. Therefore only even-parity states \k c ,n) (S6) are needed to calculate S aa {kb = 0,w). 

The calculation of dynamical structure factors of pyrochlore quantum spin ice is carried out in the same manner. 
The corresponding formula is 

- -Im5 QQ (k, W ) = --Im(k,^ a | 1 |k,^ a ) (S22) 

7T 7T U! + l0 + — H c g 

where H q h is the effective Hamiltonian for a single string in pyrochlore quantum spin ice. When the momentum 
transfer k || c, only the sum S aa + S bb contributes to the total scattering intensity due to the transverse polarization 
factor of neutrons. The relevant initial states are 

\4>a) = i(e ik - R »|k,b ) +e ik ' Rl |k,b 1 ) -e lk ' R2 |k,b 2 ) - e lk ' R3 |k, b 3 » (S23a) 

|^ 6 ) = i(e l ( k - R « +7r / 3 )|k,b ) -e^'^-^l^bi) +e i(k ' R2 - 7r/3) |k,b 2 ) - e l(k - R3+7r/3) |k,b 3 )) (S23b) 

for 5 aa and S bb respectively. Here Ri are spatial coordinates of the four inequivalent sites in pyrochlore lattice. The 

phase factors are due to the mismatch between the local spin frames (x, y, z) and global spin frame (a, 6, c). 

I 



